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ABSTRACT 

The constitutive properties of a cross-linked single-walled carbon nanotube 
material are predicted with a multi-scale model. The material is modeled as a 
transversely isotropic solid using concepts from equivalent-continuum modeling. The 
elastic constants are detennined using molecular dynamics simulation. Some 
parameters of the molecular force field are detennined specifically for the cross-linker 
from ab initio calculations. A demonstration of how the cross-linked nanotubes may 
affect the properties of a nanotube/polyimide composite is included using a 
micromechanical analysis. 

INTRODUCTION 

Carbon nanotube composites have the potential to become light-weight, high- 
strength alternatives to conventional composites for aerospace applications. To 
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facilitate the development of these materials, validated multi-scale models are essential 
to identifying the relationships between the molecular structure and the mechanical 
behavior of the material [1,2]. At NASA Langley Research Center, a class of cross- 
linked nanotube materials is being developed through a collaboration of synthesis, 
modeling, and characterization capabilities. 

Cross-linking between nanotubes has the potential for increasing the mechanical 
perfonnance of the material in nanotube based composites. Cross-linking may affect 
load transfer between nanotubes and dispersion of the nanotubes. To date, work has 
focused on the development of a short, rigid linker, and constitutive behavior of the 
material has been derived with a multi-scale model [3,4]. With this type of linker, 
some separation of the nanotubes, as well as some improvement in the transverse 
shear modulus over that of the non-cross-linked nanotube bundle was observed. 

The objective of the present work is to use multi-scale modeling to extend 
previous work to investigate the effect of a more complex nanotube cross-linker on the 
overall constitutive properties of a nanotube/polymer composite. The modeling 
procedure begins with the molecular modeling of the cross-linked nanotubes using ab 
initio calculations and molecular dynamics (MD) simulations. The elastic constants of 
the cross-linked nanotube material are subsequently detennined using an equivalent- 
continuum modeling technique. The elastic constants are then used in micromechanics 
to show the perfonnance of a cross! inked nanotube/polyimide composite material for a 
range volume fractions. 


MATERIAL SYSTEM 

A model is developed for a cross-linked single-walled carbon nanotube material. 
The cross-linking material is a symmetrical di-linker which connects two nanotubes 
(Figure 1). It contains two stiff (chemically rigid) segments near the nanotubes and a 
flexible set of ether linkages in the middle segment. The atomistic structure of the 
cross-linking connections to the nanotubes is shown in Figure 2(a), and the molecular 
structure of the thermally-equilibrated system is shown in Figure 2(b). The nanotubes 



are (10,10) single- walled achiral nanotubes with a radius of 6.78 A. The atomistic 
structure includes 9 nanotubes which are 54 A long. 
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Figure 1. Variable stiffness cross-linker. 



Figure 2: Molecular structure of cross-linked nanotube material and principal 

directions, (a) Initial structure and (b) Equilibrated structure 

MULTISCALE METHOD 

To predict the overall mechanical properties of a cross-linked nano tube/polymer 
composite material, a hierarchical multi-scale modeling procedure has been employed. 
First, the details of the molecular structure of the cross-linking agent are analyzed with 
ab initio calculations. This step is important because it provides the basis for an 
appropriate chemical structure. It generates two pieces of infonnation: the molecular 
geometry of the cross-linker, and an approximation to its electronic structure. The 


molecular geometry includes chemical bond lengths and angles between atoms. The 
electronic structure is approximated as a set of partial point charges located at the 
atomic position. These electronic charges also contribute directly into the strain 
energies used in higher levels of the model. For the cross-linking agent used in this 
work, the inclusion of electronic charge is especially necessary because 
electronegative atoms, such as nitrogen, oxygen, and fluorine, are present in the 
chemical structure (Figure 1). 

The second step of the modeling is to derive a force field to be used in classical 
simulations. This step moves beyond the molecular structure to consider which 
molecular forces will describe the molecular interactions, and ultimately what 
molecular attributes will, therefore, contribute to the strain energy. The force field 
includes analytical functional fonns for the potential energy and forces which 
detennine how the bonds stretch, the angles bend, and what torsional motion occurs 
within the molecule. It also includes tenns which describe the non-chemically bonded 
molecular interactions such as approximations to the van der Waals forces and the 
electronic forces. The van der Waals interactions are approximated by a Lennard- 
Jones potential function. The electronic charges interact via Coulomb’s law. 

Third, molecular dynamics simulations utilize the force field to detennine the 
overall equilibrium molecular structure and the strain energy for the equivalent- 
continuum modeling. The atomistic description is necessary to detennine the effects 
of the molecular structure, and especially, the effects of the chemical cross-linking on 
the strain energies. With no chemical linkers between nanotubes, the system is a 
nanotube bundle with van der Waals forces between the nanotubes. The strain 
energies are detennined from the application of a displacement field to the molecular 
system. When a displacement field is applied to the molecular structure of the cross- 
linked nanotubes, there is a change in the potential energy of the molecular system 
which is equal to the strain energy in the equivalent continuum. For the cross-linked 
nanotubes, the equivalent continuum has the properties of a transversely isotropic 
solid. The source of the contributions to the strain energies can be broken down into 
the molecular interactions in the molecular dynamics simulations. Every tenn in the 
force field contributes to the potential energy, and when the system is under strain, the 
potential energy changes accordingly. 



Fourth, the constitutive properties of an equivalent continuum are detennined from 
strain energies computed by the molecular dynamics simulations. It is assumed that 
the material is transversely isotropic with its overall mechanical behavior represented 
with five independent elastic constants. Finally, and by way of application to 
composites, the constitutive properties are used to predict the mechanical properties of 
a polymer/nanotube composite using micromechanics to estimate the effect of 
embedding the cross-linked nanotubes as an effective fiber in the polymer matrix. 

Quantum Chemistry 

The cross-linker can be thought of as a symmetric molecule composed of two 
segments as shown in Figure 1. To reduce the size of the quantum ab initio 
calculations, charges were fit for each of the two segments separately, assuming 
hydrogen tennination. The molecular geometry of each segment was optimized using 
a HF/STO-3G calculation [5]. This calculation minimizes the energy of the molecule 
in the Hartree-Fock (HF) approximation to the Schroedinger equation. The Hartree- 
Fock approximation does not include electron correlation tenns, and it assumes only 
some instantaneous electron correlation. The STO-3G refers to the basis set which is 
the selection of atomic orbitals and their analytic fonn. It includes the s- and p- 
orbitals of the first row atoms of the periodic table and hydrogen. To fit the charges at 
this geometry, a larger basis set, the 6-31 G*, was used. This basis set is extended to 
include some d-type polarization functions for non-hydrogen atoms [5]. The charges 
were fit using the ESP method [5]. 

Having obtained charges for each atom, some modifications were made to 
combine the two segments to fonn the cross-linker. First, the charge on C-atom on the 
cross-linker attached to the nanotube was neutralized. This change was to avoid an 
unnecessarily polar bond (large charge separation) as the attached hydrogen was 
replaced by a nanotube atom which was assigned to be neutral. Second, at the end of 
the propyl groups, it was observed that the aliphatic carbon furthest away from the 
C=0 bond had a large negative charge on it of -0.4 e. To fonn the hexyl group, this 
carbon, the adjacent one, and their associated hydrogens were neutralized. Minor 



adjustments were made in the rest of the linker to accommodate these changes. All the 
ab initio calculations were performed with NWChem [6]. 

Force Field 

The nanotube force field is based on the AMBER force field [7]. It includes 
energy contributions from the bond stretches, angle bends, torsions, non-bonded 
interactions between the 1-4 torsional interactions and non-bonded interactions with 
all other atoms in the system. The bonds are represented by hannonic stretches by 

U=\k b (r-r,f (1) 

where U is the potential energy, kb is the force constant for the bond stretch, and ro is 
the bond length. These values are k/,=3924 kj/mol [7] and ry= 1 .42 A. Similarly, the 
angle bends are represented by hannonic valence angle potentials of the fonn 

U = H a (O-O 0 ) 2 (2) 

where k a is the force constant for the angle bend and 6o is the equilibrium angle, 
parameterized as k fl =570 kJ/mol and 0 ( j= 120° [7]. The torsions are represented as 

U = A[ 1 + cos(n y/ - d )] (3) 

where the torsional constant A= 15.2 kj/mol, the phase shift <7=180°, the integer n = 2 
and the torsional angle is y/ [7]. Non-bonded interactions between nanotubes are 
represented with the Lennard -Jones potential. The parameters for this potential are 
given in Table I. Non-bonded interactions between 1-4 torsional partners are also 
included, but are scaled by 50%. 

Similarly, the cross-linker and its chemical bonds to the nanotube are also 
represented with the AMBER force field [7], Most of the parameters are taken 
directly from Ref. 7. In addition to the potential fonns used for the nanotube, this 



force field also includes improper torsions with the phenyl rings of the cross-linker. 
These interactions are represented with Eq. (3), and the parameters from Ref. 7. The 
cross-linker to nanotube bond was represented as by an aromatic to aliphatic C-C bond 
with ro=1.53 A. Some modifications to the force field parameterization were made. 
First, the Lennard-Jones interactions were simplified to include the list in Table I. 
This change was for convenience in setting up the force field manually in DL-POLY. 
Second, the partial point charges, derived for specifically for the cross-linker, from the 
ab initio calculations (described above) were used a classical approximation to the 
distribution of electronic charge within the molecule. While it is typical for AMBER- 
type force fields such as the one in Ref 7 to include partial point charges at each atom 
as a classical approximation to the electronic structure of the molecule, charges are 
usually specific to the local chemical arrangement of molecule, and are, therefore, best 
derived for the particular molecular structure of interest. The charges contribute to the 
system potential energy, and therefore, the strain energy, via Coulomb’s law. Third, 
the bond lengths ro and angles do for the cross-linker where chosen from the optimized 
geometry obtained from the ab initio calculations used in the charge-fitting procedure. 


TABLE I. Lennard-Jones Parameters 


Atom Type 

Lennard-Jones £>./ (k.l/mol) 

Lennard-Jones cr L/ (A) 

C 

0.3598 

3.40 

H 

0.0628 

2.60 

H (from N-H) 

0.0657 

1.07 

0 

0.7950 

2.98 

N 

0.7113 

3.25 

F 

0.2552 

3.12 


Molecular Dynamics Simulations 

The molecular dynamics simulations used to apply the displacement fields are 
started with an equilibrated structure like that in Figure 2(b). In these simulations, 
Newton’s equations of motion are integrated to move atoms according to the 
intennolecular forces. The molecular structures evolve in increments of time steps 




(0.5 fs). A strain level of 0.25 % is applied incrementally every 10,000 (5 ps) steps of 
the simulation (0.5 fs per time step). The strain is applied in the fonn of displacement 
to all the vectors describing the box shape, as well as to all the atomic coordinate 
vectors. The molecular structure was strained in both positive and negative directions 
making a total of 10 simulations to obtain five independent elastic constants. Each 
constant was then calculated by a parabolic fit to average potential energy over the last 
3,000 steps per strain level versus strain. 

To obtain the structure of Figure 2(b), molecular dynamics simulations were 
started with the initial structure in Figure 2(a). After a few steps to begin the 
equilibration process, the size of the molecular structure was gradually reduced in the 
transverse directions (2 and 3 in Figure 2). This compression was done by applying 
strain at a rate of 0.25% per 250 steps. These simulations were perfonned at a 
constant temperature of 300 K in the microcanonical ensemble (constant NVE), only 
deviating to apply the displacement field which forces a change in the system volume. 
A second simulation was used to refine the structure before applying the displacement 
fields. In this simulation, the structure was further equilibrated in an isobaric- 
isothennal (NPT) ensemble, where all three dimensions were adjustable. The 
simulations were carried out with DL-POLY [8], using periodic boundary conditions. 
An Ewald summation was used to remove the effects of the long-range forces on the 
periodicity of the structure, which result from the partial charges at each atom. This 
step doubled the computational time. 

Elastic Constants 

The material system is modeled as a transversely isotropic solid. The constitutive 


relation is 




~ Qi^n 

+ C X2 S 22 

+ ^12^33 

^22 

= C 12 £ n 

+ C 2 2 s 22 

+ C 23 £ 33 

(J 33 

= C £ 

'-'32 c ll 

+ C £ 

T <^23 c 22 

+ C £ 

-T ^22 c 33 

^23 

= ^23/23 



^12 

II 

P 

to 

to 



°"l3 

= G uYn 




where ay Cy, and Sj are components of the stress, stiffness, and strain tensors 
respectively, yy the shear strain, and i,j= 1-3 [9], The Gn and G 23 are the longitudinal 
and transverse shear moduli, respectively. Isotropy exists in the 2-3 plane, therefore, 

^12 = HQl +2^12 +^ 22 ) 

K 23 =\(C 22 +C 23 ) . (5) 

g 23 =Hc 22 -c 23 ) 

where Kn is the longitudinal plane strain bulk modulus and K 2 j is the transverse plane 
strain bulk modulus. The five independent constants Cu, K 12 , K23, Gn, and G23 are 
computed by employing displacement fields in molecular dynamics simulation. 
These displacement fields u(x) are summarized in Table II, along with the respective 
strain energy A over the system volume V. Unless the amounts of strain s and shear 
strain y are specified, the strains are zero. From the five independent constants, Cu, 
Cj2, C22, C23, G23, and G’/ jof Eq. (4) can be detennined. 



TABLE II. STRAIN ENERGY AND DISPLACEMENT FIELDS 


Elastic 

Constant 

Strain Energy 
Equivalent 

Strain 

Displacement 

Fields 

C n 

A = }VC n s 2 

£l,=S 

li(x I )=£Xi 

u(x 2 )=0 

u(x 3 )=0 

K l2 

A = 2 VK n s 2 

£ii = £22 = £ 

u(x 3 )=£Xi 

u(Xt)=£X 2 

u(x 3 )=0 

k 23 

A = 2 VK 22 s 2 


u(x])=0 

u(x 2 )=£X 2 

u(x 3 )=sx 3 

G n 

A = }VG i2 / 2 

II 

u(xi)=yx 7 
u(x 2 )=yxi 
u(x 3 )=0 

g 23 

a = \vg 22 y 2 

y. 23=r 

u(xi)=0 

u(x 2 )—yx 3 

u(x 3 )=yx 2 


Micromechanics 

The constitutive behavior was predicted with a micromechanical analysis 
method. The composite material was assumed to consist of a cross-linked nanotube in 
a polyimide matrix. It was assumed in the micromechanical analysis that perfect 
bonding exists between the nanotube/polymer effective fibers and the surrounding 
polyimide matrix. 

The micromechanics-based Mori-Tanaka method [10,1 1] was used to predict the 
elastic mechanical properties of the composite material. For this method, the overall 
elastic-stiffness tensor of the composite containing transverse-isotropic effective fibers 
embedded in an isotropic matrix material is 

C = f„C" + c f (C / T / ))(c ffl I + c f (t / )}' 1 (6) 




where c/ and c m are the effective fiber and matrix volume fractions, respectively, C 1 
and C m are the stiffness tensors of the effective fiber and matrix, respectively, I is the 
identity tensor, the angle-brackets indicate an effective-fiber orientation average, and 
T^is the dilute strain-concentration tensor of the effective fibers, and is given by 


T' = 


i+s / (c m )“ 1 (c / -c m ) 


(V) 


where S/ is the Eshelby tensor [12]. For three-dimensional randomly oriented 
effective fibers, the orientation average of a tensor, A, is 



( 8 ) 


where ij,m,n = 1,2,3; the indicial summation convention is used; 5,y is the Kronecker 
delta; and 
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Therefore, from Eqs. (8) and (9), (A) is isotropic. 

For the cross-linked, nanotube fiber/polyimide composite considered in the 
present study, the elastic stiffness components, volume fraction, length, and 
orientation of the effective fiber were used for the inclusion properties in Eq. (6). The 
effective fibers were assumed to have a spheroidal geometry for the Eshelby tensor. 
While the nanotube and effective-fiber lengths are equivalent, the nanotube volume 
fraction was detennined to be 62.5% of the effective-fiber volume fraction. This value 
was calculated by assuming the nanotube volume is defined as the total space 
occupied by the nano tube, including half of the van der Waals separation between the 



nanotube and polymer. It was assumed that the polyimide matrix was isotropic with 
properties equal to those of LaRC-SI [13], with a Young’s modulus and Poisson’s 
ratio of 3.8 GPa and 0.4, respectively. The overall composite stiffness was calculated 
for randomly-oriented effective-fibers with aspect ratios (length/diameter) of 1000 and 
volume fractions up to 40%. 


RESULTS AND DISCUSSION 

The elastic constants for the cross-linked nanotubes are included in Table 1. The 
lowest values, K23, Gn, and G23 have the highest degree of uncertainty. The reason for 
the uncertainty is that these values are the most affected by the thennal noise in the 
simulation that arises when the potential energy change resulting from the 
displacement field is small compared to the potential energy change resulting from the 
thennal fluctuations. The uncertainty in the simulation results for C11 and K12 is much 
smaller, and can be estimated from standard error of the fit. For K12, a standard 
deviation of +/- 6 GPa is estimated from the variation in the parabolic fit of strain 
energy versus strain. With this degree of uncertainty in mind, it is assumed that the 
five independently determined constants in Table III do not give a positive definite 
matrix (C12—I6 GPa). To within 95 % confidence, K12 is within a range of 68-92 
GPa. An arbitrary choice of Ki2=88.6 GPa is chosen to obtain a positive definite 
matrix. This value is reasonable within the simulation uncertainty. With this value, 
the six constants in Eq. (1) become C11 = 350 GPa, C22 =3.5 GPa, Ci2= 1.1 GPa, C23 = 
1.9 GPa, G12 = 3.7 GPa, and G 23 = 0.8 GPa. 

For comparison, the elastic constants of a pure nano tube bundle (without cross- 
linker) with the same nanotube force field [3] are examined. The current C11 value is 
61 % of the Ci 1 value of the nanotube bundle which is close to the proportionality 
expected from the reduction in nanotube volume fraction (62.5% in the cross-linked 
nanotube fiber). The value of G23 also compares favorably. The current value of G12 
is an order of magnitude lower, possibly because the flexible cross-li nk er is able to 
absorb some strain energy. The current value of K23 is also an order of magnitude 
lower, most likely because of the reduction in van der Waals attractions from the 



nanotube separation in the cross-linked material. This effect was observed with a 
shorter cross-linker [3]. The most unusual difference is in the K 12 value. This increase 
may simply be the result of packing a large amount of cross-linker between the 
nanotubes, therefore, contributing additional resistance to defonnation in both 
directions. 


TABLE III. ELASTIC CONSTANTS 


Elastic 

Cross-Linked Nano tube 

Nanotube Bundle Without 

Constants 

Fiber (GPa) 

Cross-Links (GPa) [3] 

Cn 

350 

570 

Kj2 

80 

11 

K 23 

2.7 

24 

G 12 

3.7 

35 

G 23 

0.84 

1 


The Young’s and shear moduli of the embedded, cross-linked nanotube effective 
fiber in the polyimide LaRC-SI are plotted in Figure 3 as a function of nanotube 
volume fraction for randomly orientated libers. Both the Young’s and shear modulus 
show monotonic increases with increasing nanotube volume fraction. 
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Figure 3. Elastic constants of the cross-linked nanotube/LaRC-SI Composite. 


CONCLUSIONS 

The constitutive properties of single-walled carbon nanotubes cross-linked with a 
variable stiffness tethers were detennined with a multi-scale modeling method. The 
method used molecular dynamics simulation and equivalent continuum modeling to 
detennine the elastic constants. Input for the force field used in the molecular 
dynamics simulations was customized for the cross-linker with ab initio calculations. 
It is assumed that the material was transversely isotropic, and therefore, five 
independent constants were determined. A comparison of the constants with similarly 
derived ones for a nanotube bundle without any cross-linker demonstrated that the 
mechanical properties are affected by the separation of the nanotubes in the bundles, 
and by the cross-linker having multiple degrees of freedom giving it the capability of 
absorbing some of the strain energy. The effect of the cross-linked nanotubes in 
increasing the Young’s and shear moduli of a composite with increasing nanotube 




volume fraction is demonstrated by modeling the nanotubes as effective fibers in a 
polyimide with a micromechanical analysis. 
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